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Abstract. In a previous contribution (H.J. Stockmann, J. Phys. A35, 5165 (2002)), 
the density of states was calculated for a billiard with randomly distributed delta-like 
scatterers, doubly averaged over the positions of the impurities and the billiard shape. 
This result is now extended to the fc-point correlation function. Using supersymmetric 
methods, we show that the correlations in the bulk are always identical to those of the 
Gaussian Unitary Ensemble (GUE) of random matrices. In passing from the band 
centre to the tail states, The density of states is depleted considerably and the two- 
point correlation function shows a gradual change from the GUE behaviour to that 
found for completely uncorrelated eigenvalues. This can be viewed as similar to a 
mobility edge. 



PACS numbers: 05.45.Mt, 03.65.Nk, 05.30.-d 

Submitted to: J. Phys. A: Math. Gen. 
1. Introduction 

The theory of random matrices provides a schematic, but powerful statistical model for 
a wide class of spectral problems in complex systems, for reviews see Refs. El E] • In 
particular, there is overwhelming evidence for the fact that the spectral fluctuations 
of a quantum system whose classical counterpart is fully chaotic are described by 
the Gaussian ensemble of random matrices, i.e. by the Gaussian Unitary Ensemble 
(GUE) in the absence of time reversal invariance and by the Gaussian Orthogonal 
Ensemble (GOE) if time reversal invariance holds and the spectrum is free of Kramers 
degeneracies j3]. On the other hand, the fluctuation properties for quantum systems 
whose classical counterparts are regular ought to be different, and often of the Poisson 
type. Many systems show mixed fluctuation properties and transitions from regular to 
chaotic behaviour. 

Quantum billiards are ideal systems for the study of spectral fluctuation properties. 
Billiards are said to be ballistic because the classical dynamics and the quantum spectra 
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are exclusively determined by the shape of the boundary Whereas such ballistic systems 
are well understood, the situation is less clear for disordered systems. In particular, there 
are many open questions concerning the localization-delocalization transition if disorder 
is varied. From the one-parameter scaling hypothesis [5] is is generally accepted that 
in one- and two-dimensional systems all states are localized, but analytic proofs exist 
only for one-dimensional systems (see reference [S] for a review) . There are a number of 
work using supersymmetric techniques, where the change of the wave function amplitude 
statistic is studied with the reciprocal conductance as a perturbation parameter [Zj, but 
up to now there is no closed theory covering the full range from localized to delocalized 
wave functions. On the other hand there are microwave experiments showing a clear 
localization-delocalization transition with frequency [HE]- 

This was the motivation of a previous publication ^Qj, hitherto denoted Ref. I, to 
tackle the problem by an alternative approach. Instead of the usually applied non-linear 
cr-model the more explicit system of a billiard with randomly distributed scatterers was 
studied. This approach generalized a model introduced by Bogomolny et al. jTTj. The 
average over disorder was achieved with help of a trick using the conjecture that a typical 
wavefunction can be viewed as a random superposition of plane waves ^2] • Thereby, no 
supersymmetric field variables are needed which Efetov used to construct his non-linear 
cr-model [T^j. It avoids as well the complications of diagrammatic expansions of Green 
functions and summations of ladder diagrams 6j. It was already conjectured in Ref. I 
that there should be a localization-delocalization transition with increasing number 
of scatterers. In the present work further arguments are given that for a sufficiently 
large number of scatterers there is indeed a mobility edge, separating the band from 
the tail states, where such a transition takes place. These effect is accompanied by a 
considerable depletion of the density of states. There is a fundamental difference to the 
cr-model which will be discussed. 

The article is organized as follows. In Section 121 the main results of Ref. I are 
recapitulated and the /c-point correlation function is calculated generalizing a method 
developed in reference [15J. In SectionElthe results are specialized to the strong coupling 
limit, and it is shown that everywhere within the band random-matrix results are 
recovered. In Section H] the behaviour of the /c-point correlation close to the band 
edge is studied. The two-point correlation function in particular shows a transition 
from GUE behaviour to that of completely uncorrelated eigenvalues suggesting that 
there is indeed a mobility edge. 

2. The Model and its Supersymmetric Evaluation 

We setup the model in Sec. 12.11 and map it onto superspace in Sec. 12.21 The kernel 
determining all correlation functions is calculated exactly in Sec. 12.31 The density of 
states is worked out in Sec. 12.41 in the strong coupling limit. 
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2.1. Setup of the Model 

In Ref. I the density of states was calculated for a billiard with randomly distributed 
scatterers, averaged of the the positions of the scatterer. The system was described by 
the Hamiltonian 

H = H + V, (1) 

where Hq is the operator of kinetic energy, and V is the scattering potential. Assuming 
L point-like scatterers at positions ri, we have 

L 

V(r) =4 7 rA^5(f-f i ) . (2) 

1=0 

Using standard supersymmetric techniques, the density of states was expressed as the 
derivative 



(3) 

j=o 



of the generating function 

Z(E 1 ,E 2 )= (4) 
J dfx] exp \ i^2[(E 1+ 6 a p - (#„)«/?) x*aXf> + (E 2+ 5 a/3 - (H o ) a0 ) j M L , 

N 

with the volume element d[x] = T7 dx*dx a d£*d£a,. Here, the quantity M is given by 

a=l 

M=(ex V ^-4Ki\J2r a (r)Mr)(x>p + € a ^ , (5) 

where the brackets denote the average over the scatterer position. To perform the 
average, in Ref. I a trick was applied by replacing the average over the positions by 
an integral over the wave function amplitudes if) at the positions of the scatterers with 
the amplitude probability density p(ip) as a weight function. For the latter a Gaussian 
distribution was taken typically for chaotic billiards [T2llllj. In the next step, a second 
average was performed by replacing the billiard spectrum with that of a random matrix 
from the GUE of rank N. 

As a result a simple analytic expression was obtained for the density of states. 
For L > N a qualitative change in the density of states was observed suggesting 
a localization-delocalization transition. In the following the results of Ref. I will be 
generalized to the calculation of the /c-point correlation function, and further evidence 
will be presented of the existence of localized states and a certain type of mobility edge 
within the present model. 

2.2. Supersymmetric Matrix Model 

To compute the fc-level correlation functions of k energies E p , p = 1, . . . , k, we combine 
and extend the procedures outlined in Refs. I and [T5] . We construct the functions 
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Rk(Ei, . . . , Ef.) obtained from averaging over the product of k Green functions, including 
their real parts. For example, the density defined in Eq. (J3J) is the imaginary part of 
Ri(E). The correlation functions Rf.(Ex, . . . ,E k ) for the imaginary parts only can be 
calculated as proper linear combinations from the functions R k (E\, . . . , E k ). The latter 
are given as the derivatives 



of the generating function 

( N 



(6) 

j=o 



(Z k (E + J)) = 2 fc ( fc ~ 1) J d[S}exp (^—TtS 2 + TtS(E + J) 

Det N S 



(7) 



Det L (l 2fc + AS) 

with respect to k source variables J p , p = l,...,k. Energies and sources 
variables are ordered in diagonal matrices E = diag (Ei, E-^, . . . , E k , E k ) and J = 
diag (— Ji, +Ji, . . . , — Jfc, + Jfe). The generating function is the straightforward 
extension of the generating function used in Ref. I to arbitrary k. To keep with the 
notation in Ref. I, we introduced the 2k x 2k Hermitean supermatrix S which can 
be mapped onto the supermatrix a used in Ref. ^Hj by exchanging its bosonic and 
fermionic eigenvalues. Moreover, we use the symbols Tr and Det to indicate supertrace 
and superdeterminant. 

As in Ref. jT3] , the supersymmetric extension of the Itzykson-Zuber integral can be 
employed to reduce the generating function to an integral over the fermionic and bosonic 
eigenvalues is P 2, p = 1, . . . , k and s p x, p = 1, . . . , k, respectively. This is so because the 
term coupling S and E + J is the only one in the integrand which is not invariant 
under an unitary transformation of S. Again, as in Ref. [ISj, an easy evaluation of the 
derivatives with respect to the source variables is possible and we arrive at 

1 f,u^„J^ 



R k (E h ...,E k ) = j—— J d[s}B k (s)exp (_Trs 2 + Tr S £) 



Det N s 



Det L (l 2fc + Xs) 



(8) 



We collect the eigenvalues in the diagonal matrix s = diag (isu, . . . , is k 2, sn, . . . , Sfci). 
We notice that the bosonic eigenvalues carry a small imaginary increment, s pl = s p i + ir], 
where it is necessary. It is send to zero at an appropriate point of the calculation. 
Keeping this in mind, we can integrate all eigenvalues over the entire real axis. This is 
equivalent to the choice of the integration contour in Ref. I. The function Bk(s) 

B k (s) = det ^— (9) 

S p i ' lS q2j ptq=1 ^ 

in Eq. (JJJ) is the square root of the Jacobian which is due to the change of variables from 
the Cartesean coordinates in S to eigenvalues s and angles. It is a determinant which 
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couples one bosonic and one fermionic eigenvalue in each of its elements. Expanding 
the determinant, 

g*(«) = Se(7r)]I - _) , (10) 

n 1 = 1 °pl 6t, qn(l) 

where the sum is over all permutations tt, and e(7r) = ±1 for even, and odd permutations, 
respectively, the integrations in equation (|SJ) factorize into products of double integrals, 
each over one bosonic and one fermionic variable. The result can again be written in 
terms of a determinant 



R k {E 1 ,...,E k ) = det C NL {E p ,E q ) (11) 

L ip,q=l,...,k 

with a kernel given by 

a , ^ ^ x 1 f + °° r +oc dsids 2 ( N , 9 9 , 




+ts ^- SiE '){TTMr^ I -I ■ (12) 

We suppress the indices p and q in the integration variables. Thus, the correlation 
functions have a determinant structure which is a immediate consequence of the 
determinant (JHJ). In full analogy to Ref. [T^], we obtain the correlation functions 
Rk(E\, . . . ,E}~) by replacing l/(sf) N in Eq. ()12|) with its imaginary part Im 1/(5+)^. 
As in Ref. I, we rescale the energies and the strength parameter according to 



71 JN/2 
e v = E v and a = - — - — . (13) 

On this scale, the correlation functions are given by 

R k {e 1} ..., e fc ) = det [C NL {e p , e q )} pq=K k , (14) 



where the kernel now reads 



1 f+°° f+°° dsids 2 



V 9/ 2 

7T J— oo J —oo S\ — 



exp f — (sj + S2) + 2is 2 £g - 2si£ p 



^rim^. (15) 



,. CX -\- IS2 , 

Due to the rescaling ()13j) . we obtain the kernel for the GUE correlation functions exactly 
in the form given in Ref. jT3], if we consider the limit A — > 0, i.e. a — > 00, or, equivalently, 
L = 0. We notice that the scaling factor in Eq. (|T3~j) is precisely the GUE mean level 
spacing ix/ylN in the center of the semicircle. 

In this derivation, we have omitted a Efetov-Wegner or Rothstein contribution 
HEllIZllinj which adds to the real part of Cnl(E p , E q ). The functions Cnl(E p , E q ), the 
main objects of our interest, are not affected. 
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2. 3. Exact Computation of the Kernel 



Extending the methods of Ref . [TH] , the kernel can be evaluated exactly for all values of 
N, L and a. We define the functions 



-2i) N a L /•+«> / . . , 2 x s% 
-= — / ds a exp(-(s 2 -*e) 7 — 

y 7T J-oo v ' (a + IS 2 ) 



vnl{£) = r exp (e 



exp (e 2 ) ^-j; exp + ef) (a + Sl ) L 



(16) 



/+°° / \ 1 

dsi exp (-(si + e) 2 J (a + si) L Im 

= — ^ W 

which reduce to the Hermite polynomials Hn{s) for L = or, equivalently, for a — > oo. 
In Appendix A[ some properties of these functions are compiled. We now express 
lml/(sf) N in Eq. (fTB} as d N ~ 1 S(si)/dsi _1 and integrate by parts until the (N — 1)- 
fold derivative with respect to Si acts on all Sj dependent terms in the integrand. After 
applying Leibniz' rule for multiple derivatives of products, we can insert the second form 
of the function vnl(e) into Eq. ([T5|). The S2 integration then yields just the function 
unl(z) and we arrive at 

1 N ~ x 1 

CVl(£ p , e g ) = -f= exp (-ej) ^ t^-t^l^pKl^) • (17) 

V ^ n=0 A U - 

Thus, we have expressed the kernel and all correlations in terms of the functions vnl(^p) 
and UNL^Eq). Formula (fTTj) is a generalization of the corresponding expression for the 
GUE. We mention in passing that one also derives 

C NL (e p ,e q ) = 2N - 1{N _ 1)l ^ 

/ exp(-(e q + t) 2 ^u NL (e q + t)v {N _ 1)L (e p + t)dt , (18) 

•/ 

which again generalizes the corresponding expression for the GUE in Ref. |15| . The 
result (fTSj) involves only the orders A" and A" — 1 of the functions vnl(e p ) and UNL(s q ), 
which are not even orthogonal. 

2-4- A Christoffel-Darboux Formula for Cnl(e p , e q ) 

For L = or, alternatively, a — > oo, the sum on the right hand side of expression (|17|) 
can be performed with the result 

1 

2 Jv-i(AT~ i)!0f^V -i) 

U(N-l)o(£q)VNo(£p) — UN ( £ q) V (N-l)o(£p) 



C m (e p , £ q ) = oN _ UM — -tt^ exp (-E 2 



E p Eg 



(19) 



This is the well-known Christoffel-Darboux formula for the Hermite polynomials. This 
expression is now generalized to arbitrary values of L. To this end we multiply both 
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sides of equation ()15|) by e p — e q and obtain in a sequence of elementary steps, including 
one integration by parts, 



(Ep — £q)CNL(£p, Eq) 



-Im 



7T 



00 /-+ 00 dsids 2 



00 S\ — IS2 



exp (—(si + si) + 2is 2 e q - 2siE 1 



a + s 



a + is2 



is 2 

A, 



N 



— Im 

7^ 



00 /-+ 00 dsids 2 



X 



d 



00 Si — 2S 2 

1 d 



2 I 9s 1 i ds 2 



x exp ( — (si + s 
1 



exp (2is 2 £ q - 2s 1 e p ) 

N 



a + si 



-Im 



7T 



a + is 2 
00 /-+ 00 dsids 2 



is 2 

A, 



X 



<9 



00 S\ — IS 2 

1 d 



exp (2is 2 £ q — 2si£ ? 



2 \<9si i <9s 2 



1 u \ f . 2 



OJ + Si 

a + is 2 



2S 2 



AT 



1 /.+00 r+00 dsids 2 
— Im / / — 

7T Z J-00 Si — 2S 2 

1 / L 
x - [ -2si + 2is 2 + 



L N N 

a + Si a + is 2 is 2 si 



x 



exp (—(si + si) + 22S2£ g - 2si£ ? 



a + si 
a + is 2 



«s 2 



N 



rim 

1 



00 r+00 



-2 - 



ds\ds 2 
L 



+ 



N 



(a + si)(a + is 2 ) is 2 sj 



x 



exp (—(si + si) + 2is 2 e q - 2si£j 



a + si 



- , ^ ■ (20) 
.a + zs 2 / \Si / 

The inconvenient denominator coupling the si and the s 2 integrations has disappeared 
with the consequence that all integrals can be expressed in terms of the unl(^) and 
vnl(e)- A formula of the Christoffel-Darboux type obtains 

1 



£q)CNL(£p, £q) 



2 N (N-1)\^ 



exp 1 —e 



~ U NL (£ q )V( N ^ L (£ p ) + U( N „ 1)L (£ q )v NL (£ p ) 



^2 Un ( l + 1 ) ( £ l) V (N-l)(L-l) (e P ) 



(21) 



which is valid for all values of N, L and a. This is quite remarkable, because the 
functions u^l{e) and Vnl(e) are no orthogonal polynomials. In different context, similar 
generalizations of the Christoffel-Darboux formula have been obtained in Refs. [2D 122] 
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and by Strahov and Fyodorov ^Sj i n the calculation of correlation functions of ratios 
and products of characteristic polynomials of Hermitian random matrices. 

For L = expression ()19j) for the Hermite polynomials is recovered. Another special 
case is obtained for the strong coupling limit. Here the Gauss functions in the integral 
(fTKj) may be replaced by one with the consequence that the first term in the brackets 
on the right hand side of equation (j21j) is missing. Furthermore in this limit the unl(s) 
and Vnl(e) can be expressed in terms of generalized Laguerre polynomials, 

u NL (e) = v ^exp(£ 2 )(-l) iV (2a)( iV+1 ) 

x exp(-^-^Lr N ^(z) (22) 

v NL (e) = (-l) N N\a- N L%- N) (z) , (23) 
where z = lea = E/X as in Ref. I. Collecting the results we obtain from equation (J2~T|) 



(s p - e q )C NL (e pj e q ) = — — — - exp 



-z n )z 



L-N 



x [-L^\z q )L^ N \z P ) + Lr N \z q )L^\z P )\ • (24) 

Comparing equations (J17j) and (J2ip we obtain the following Christoffel-Darboux relation 
for the generalized Laguerre polynomials 

JV-l 

Y,n\x N -^L^\x)L^\y) = 

n=0 

j\t\ Ln \ x )^n-i \v) ~ ^n~i \ x )^n \v) ^5) 
y-x 

This is not the Christoffel-Darboux relation for the Laguerre polynomials found in 
compilation such as Ref. |19j . but we cannot exclude that it is known in the mathematical 
literature. 



3. Density of States and Correlations via a Saddlepoint Approximation 

In Section 13.11 we work out the density of states in the strong coupling limit. The 
correlations in the bulk of the spectrum are computed for arbitrary coupling in 
Section 13.21 



3.1. Density of States in the Strong Coupling Limit 

For strong coupling A > 1 or a C 1 and L > N, the density of states was evaluated 
in Ref. I by means of a WKB approximation to leading order in L > N ^> 1. Here we 
show that this is equivalent to a saddlepoint approximation. For k = 1, we write the 
generating function (fTTj) in the form 

(Z 1 (e + J))= J d[S] exp (jC(S,e + J)) 

C(S,e + J) = TrS 2 + 2TrS(e + J) + NTrlnS - LTr\n(a + S) . (26) 
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2ae 2aA(e-e_) 



Figure 1. (a) Density of states in the strong coupling limit for L/N = 4. (b) Density 
of states in the neighbourhood of the lower band edge (for details see section ET2t . 



We use the rescaled variables ()13j) . drop the index 1 on the energy variable and write 
e shorthand for el2- In the strong coupling limit S is of the order of a as can be seen 
by applying the substitution S = aS'. The term Tr S 2 is of thus of the order of a 2 and 
may be dropped in the Lagrangean C(S, e + J). In this approximation, the saddlepoint 
equation resulting from the condition dC = at J = reads 
N L 

2e + = 0, (27) 

s a + s 

with s standing for the two scalar saddlepoints Si and is 2 o- The solutions can be 
written as 

s = — (-{2ea - (L - N)) T i^A-LN - (2ea - (L - N)) 2 ^j . (28) 

Obviously, the imaginary part is only non-zero if the energy satisfies 

S-<e<e + with s T = — (L + iV =f 2VZ~n) . (29) 

2a ^ ' 

We now expand the Lagrangean C(S, e + J) around the saddlepoints up to second order 
and integrate out the massive modes in a Gaussian fashion. One can convince oneself 
in a straightforward, but tedious, calculation that these Gaussian integrals converge as 
long as the condition (|2"TI|) holds. At the saddlepoints, the Lagragean is simply 4so</i 
and we find from Eq. © 

2 

Ri(e) = -s 

7T 

1 



-(2ea - (L- N)) + i^ALN — (2ea — (L — N)) 2 ) . (30) 



2txe 

This is the full one-point function in the strong coupling limit. The imaginary part is 
the density of states which is non-zero for e_ < e < e + . As expected, it coincides with 
the WKB approximation of Ref. I. The saddlepoint approximation yields, in addition, 
also the real part of the one-point function. As an illustration figure ^a) shows the 
density of states calculated from the imaginary part of Ri(e) for L/N = 4. 
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3.2. Correlations in the Bulk of the Spectrum 

We take advantage of a remarkable identity which connects the kernel and the generating 
function for k = 1, 

C NL {e p ,e q ) = — AZ X {B)) 



7T[£ p - E qJ 



7r(e p - e 



J d[S] exp (TtS 2 + TtSS) 
Bet N S 



(31) 



Bet L (al 2 + S) 

with £ = di&g(e p ,e q ). This identity which is easily derived with the supersymmetric 
extension of the Itzykson-Zuber integral for k — 1 allows us to work out the correlations 
on the unfolded scale by a saddlepoint approximation involving 2x2 supermatrices, 
i.e. in a simple Cartesian space. A similar procedure was employed in Ref. (201 i n the 
context of chiral random matrix ensembles. 

To begin with, we discuss the strong coupling limit of the previous section and turn 
to the general case later. We write 

.r_p 
2 



£ =e pq l 2 + D^A with 



e pq = £p £q and Dr pq = e p - e q . (32) 

Here, we introduced the metric A = diag(+l, —1) and, anticipating the steps to come, 
the local mean level spacing D = 1/R\{e pq ) which defines the unfolded scale. We use the 
form ()26|) for the generating function with e + J replaced by £ to evaluate Eq. ()31|) . In 
the strong coupling limit, we neglect the term TrS* 2 . Although a <C 1 in this limit, we 
do not make any assumption about its value in the present discussion. The saddlepoints 
are the stable points of the integrand in an asymptotic 1/N expansion. The unfolded 
correlations live on the local scale of the mean level spacing D. Thus, we have to keep 
r P q = (^p — £ q)l D fixed in the asymptotic expansion for the calculation of the unfolded 
correlations. The energy difference e p — e q itself appears in the integrand. As it is 
given by Dr pq , and as the mean level spacing D vanishes in the limit N — > oo, the 
energy difference cannot yield a contribution to the saddlepoints and we may neglect 
it when calculating them. Thus, we are left with exactly the same problem as in the 
previous section, only e is replaced by e pq . This implies that the integrals over the 
massive modes converge in the non-zero region of the spectrum and, moreover, that the 
only non-vanishing contribution to the correlations comes from the term TrS^ in the 
Lagrangean. Collecting everything, we find 

CNL{r pq ,r pq ) = lim DC NL (e p ,e q ) = exp(7rr pg ) eXP mrpq (33) 

for the kernel on the unfolded scale. We notice that the result depends on r pq = 
(e p —e q )Ri(e pq ) where Ri(e pq ) = Re R\{e pq ) is the real part of the one-point function. As 
discussed in Ref. [2T], an Efetov-Wegner or Rothstein term has to be added to Eq. (JBlIj) . 
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It affects only the real part and reads — l/irr pq . For the correlation functions involving 
the imaginary parts of the Green functions, we only need the imaginary part 

c N L{r pq , f pq ) = lmc NL {r pq , f pq ) = exp(7rf pg ) Sm 7rTp9 , (34) 

TTf'pq 

which consists of the GUE sine kernel and an exponential function depending on r pq . 
Both variables, r pq and r pq , are odd under the exchange of the indices p and q. Thus, 
the sine kernel stays unchanged, while the exponential function acquires a sign in its 
argument. This implies for the correlation function on the local scale 

X k (r 12 , ri 3 , . . . , r(fc_i)fc) = lim D k R k (e 1 , ...,e k ) 

L>N— +00 



= det 



sm^ (35) 



which is identical to the standard GUE correlations. 

The previous derivation is for the strong coupling limit. In the following, we present 

a general discussion of the correlations. We write the kernel as the convolution 

~ 1 r+°° / „\ r+00 . . 

C NL (e p ,e q ) = - dy p exp \-y p J / dy q exp \-y*J 

B NL {£p + y P ,£q + Wq) 
ft / n 1 / ,+ °° / ,+0 ° dSldS2 / 2 , . \ 

Bnl{z p , z q )= / / — exp -{s 1 + s 2 ) + %s 2 z q - s x z p 

TT* J-00 J-00 Si — ISo \ 



-00 S\ — IS2 

a + si \ L ( is 2 



N 



(36) 

Obviously, the kernel B NL (z p , z q ) is the kernel of the strong coupling limit. However, it 
emerges due to the convolution. Thus, we do not need to assume that L > N. In the 
following, we only assume that both numbers, L and iV are large. Moreover, we make 
no assumption about a. As we have seen in the previous discussion, this kernel, here 
denoted B NL (z P) z q ), leads to standard GUE correlations on the unfolded scale. One 
might argue that this does not necessarily carry over to the present case, because the 
arguments z p = e p + y p and z q = e q + iy q contain the integration variables y p and y q . 
However, as we are only interested in the fluctuations, we only need to consider the 
integration variables on this scale. Thus, we may neglect them for the determination of 
the saddlepoints. After assembling things properly, we arrive at 

c N L(r P q,r pq ) = lim DC NL (e p , e q ) 

rs — >oo 



I r+00 , v r+oo , . 

- J_ dy p exY>(-y 2 p j J_ dy q e^[-y 2 q ) 



exp [TT(r pq + (y p - iy q )R 1 (e pq 
exp (iTr(r pq + (y p - iy q )/D)) 



(37) 



n(r pq + (y p -iy q )/D) 

As we are only interested in the imaginary part, we may again ignore the Efetov- 
Wegner or Rothstein term. The imaginary part can be obtained from the difference 
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of a retarded and an advanced Green function. The two Green functions yield the 
same kernels, apart from a sign change in the argument of the exponential function 
in the numerator, exp (±in(r pq + (y p — iy q )/D)). Hence, only the difference of these 
two exponential functions, the sine function, enters. This is equivalent to taking the 
imaginary part of Eq. (|H7|) while formally ignoring the imaginary unit coming with the 
variable y q . Thus, we find 

CNL{r pq ,r pq ) = lmc NL (r pq ,r pq ) 

= -/ dy p exp(-y ) dy q exp(-y q 

exp [7r(r pq + (y p - iy q )R x {£ pq ))) 
sin (ir(r pq + (y p -iy q )/D)) 



rr(r pq + (y p -iy q )/D) 
„ ,smnr pq 

exp(7rr w ) , (3* 



Trr pq 



where the integrals over y p and y q were done as in Ref . [22] • Hence, the correlations are, 
once more, of the standard GUE type. Some comments are in order. First, it should 
be clear that the mean level spacing D in the calculation above was formally the one of 
the strong coupling limit and has thus to be smoothly adjusted when going into another 
regime. Therefore, our line of arguing is correct only if we are always in the bulk of the 
spectrum, i.e. far away from any possible edges or gaps. Second, the discussion beyond 
the strong coupling limit could also be done in a saddlepoint approximation of the full 
expression (jUlj) . This, however, leads to a most inconvenient third order saddlepoint 
equation. In the approach chosen here we avoid this an also gain the insight that the 
strong coupling limit and the general case are related via a convolution. Third, we 
emphasize that the connection (}3*Tj) between the kernel and the generating function for 
k = 1 simplifies the calculations enormously: the saddlepoints are isolated, no Goldstone 
modes occur. Furthermore, all correlations are treated at once. 



4. The band edges 

From the pioneering work of Mott and Anderson it is known that in disordered systems 
there are no sharp band edges for the density of states. There is a mobility edge instead 
separating the delocalized states in the band from the localized ones in the tails. The 
mathematical origin of the band edges is due to the fact that in dependence of some 
parameter the two solutions of the saddle point equation (j2*7| change from complex 
conjugate to real. This behaviour is generic, though the present model the band edges 
are only an artifact of the finite rank of the matrices. It therefore is worthwhile to study 
the regime of the band edges somewhat more in detail. 

After obtaining a WKB approximation for the kernel in Section 14.11 we work out 
density of states and correlations in Section 14.11 
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4-1- A WKB Approximation for C]yL(Ep,E q 



To keep the discussion simple, we again concentrate on the strong coupling limit. 
Starting point is the Christoffel-Darboux relation (ffifj) for Cnl(£ p , £ q ) holding in this 
limit. Using standard relations for the generalized Laguerre polynomials, it can be 
written in the alternative form 



(s p — e q )CNL(e p , E q ) — 2a 



(L 



~ Z q) Z q 



„L-N 



-(L-N-l)', . (L-N-l), x 
J N {ZqJ-^N \ Z P) 



- (L-N-l), 



(L-N-l)' 



(39) 



where z p / q = 2ae p / q , which is somewhat more suitable for the present purpose. Following 



Ref. I, we write N X \z) as 



\L-l)\ 



L-N 

exp(z/2)z 2 f(z) 



(40) 



where f(z) is a solution of 

f"(z) + q 2 (z)f(z) 







N + L 1 l-(L-iV-l) 2 
~^2z 4 + 4^2 



q 2 (z) may be written as 



Az 2 



(z - z-)(z - z + ) 



where 



z± = N + L ± JANL + 2(L-N)^N + L± 2\[WL . 



(41) 



(42) 



(43) 



We note that this is the same expression, which was obtained above from the saddlepoint 
approximation for the band edges in the strong coupling limit (see section 13. lj) . 
For z_ < z C z + the WKB solution of equation (}4~Tj) is given by 



where 



/(*) 



Q(z) 



cos 



Q(z)- 



71 



q(t) dt 



(44) 



(45) 



(To be concise we restrict the discussion to the neighbourhood of the lower edge, but it 
is straightforward to transfer all results to the upper edge as well). Inserting this into 
equation (jHSJ), we recover the result for Cjvi(%,£ 9 ) obtained in section IH~2l bv means of 
the saddlepoint technique. For 2<2_ the corresponding expression reads 

l / : 



f(z) 



exp [- \Q(z) 



(46) 



2]j7r\q(z)\ 

Inserting expression into (jSHj) one notices that Cnl(£ p , £ q ) vanishes within the limits 
of the WKB approximation applied. To describe this regime appropriately, one would 
have to go to the next WKB order. 
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We do not proceed further in this direction, but concentrate on the immediate 
neighbourhood of the lower edge which is not covered by equations (}4*4*j) and ([46)1 . 
Linearizing q(z) close to Z-, 



?(*) = iM^C*-*-). (47) 



4z: 

equation (jH| can be solved with the result 
/(*) = ^Ai [\(z- - z)) , 



a/4 \ I 



V y \n + l-2Vnl) 

where Ai(z) is the Airy function. With the factor A -1 / 2 the asymptotic behaviour of the 
Laguerre polynomials is reproduced correctly by equation (|40p. This can be shown by 
techniques described e.g. in chapter 9.3 of Ref. [23j. Collecting the results, we obtain 
from equation (J3~9"j) 



C NL (£ P ,£ q ) = 2aexp f^- Zq - Zp ^j 



L-N 
2 



f'(z q )f(z p )-f(z q )f(z p ) fT AT J{z p )f{z q ) 



(L-N)- 



Zp Zq 2,ZpZq 



(49) 



The second term on the right hand side vanishes for L — > oo, since z p , z q are of order 
0(N + L), and will be discarded in the following. Essentially the same approach to 
describe the behaviour of correlation functions close to the band edges was applied by 
Akemann and Fyodorov j2l] in the study of characteristic polynomials. 



4-2. The Density of States and the k-point Correlation Function 
The density of states is obtained from equation ()49|) 

p(e) = C NL (e,e) = 2a {[f'(2ae)} 2 - f (2ae) f" (2ae)} . (50) 
For the regime close to the lower band edge we obtain by inserting expression pHjl for 
/(*), 

p(e) = 2aA{[Ai'(-s)] 2 + s[Ai(-s)] 2 } , s = 2aX{e-e-). (51) 

Figure [T^b) shows a plot of the density of states in the transition regime as obtained 
from equation (|5*T|). 

From equation (|14|) the two-point correlation function results as 

Rkiex, . . . , e fc ) = (2a) fc det [c LN (e p , £q)} M=1 ^ k , (52) 

where 

CLN{£pj£q) = — • (Od) 

Zp Zq 
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(The first factor on the right hand side of equation (J49|) cancels in taking the determinant 
as is easily seen.) Just as in section we now introduce rescaled variables 

Ep + Sq _ 
c pq 2 ' PI P 9' \ / 

and a rescaled correlation function Xfc(ri 2 , ri 3 , . . . ,r(k-i)k) — D k Rk(ei, . . . , where 
D^ 1 = p(e) is the local density of states as given by equation (|5U|). Note that in 
contrast to section EOl we do not perform the limit L,N — > oo, since we are interested 
in particular in the behaviour close to the band edges. We then have 

X k (ru, r 13 , r (fc _ 1)fc ) = det [c LN (e pq , r pq )} p q=h k (55) 

with 

/ Dr Dr\ 

cln(£, r) = c LN e + — , e — . (56) 



2 ' 2 

Inserting for f(z) the expression ([48)1 . we obtain for the regime close to the lower edge 

eiy(e , r) = ^(-«^(-'-)-^-^(-«-), (57) 

r 

where s± = s ± aXDr. In the limit r — > we obtain clat(£,0) = 1 as it should be. 
This is a direct consequence of the differential equation Ai"(z) — ^Ai(^) = of the Airy 
function, and equation (jBTjl . For r — > oo CLAr(e,r) decays according to 

2 



CLiv(^,^) ~ r 1/4 exp 



~ («ADr) 2/3 



(5? 



which follows from the asymptotic behaviour of the Airy function (see also equation 
(J46|) ) . The transition between the two regimes is observed at r = s/aXD. With 
decreasing density of states p(e) — 1/D the transition point is thus approaching r = 0, 
i.e. the eigenvalues become more and more uncorrelated. This is illustrated in figure 121 
where the two-point correlation function i?2( r ) — 1 — [cla^; r )] 2 is shown for different 
values of s — 2a\(e — e_) in the neighbourhood of the lower band edge. In addition the 
GUE result is shown for comparison. We observe with decreasing s a gradual transition 
from a GUE behaviour to that expected for completely uncorrelated eigenvalues. This is 
exactly which is expected for a mobility edge: within the band the eigenvalues experience 
a quadratic level repulsion typically for the GUE, whereas in the tails the localization of 
the wave function leads to a suppression of the level repulsion (see e.g. reference [25J). 



5. Summary and Conclusion 

In Ref. I the density of states for the billiard with randomly distributed scatterers 
was calculated, doubly averaged over disorder and shape of the billiard. We mention 
in passing that the resulting model shows some formal similarities to chiral random 
matrix models [26J. This is due to the way how the average over the disorder is done. 
In the present work, the results of Ref. I are extended. We calculate the fc-point 
correlation functions exactly. The model of Ref. I generalizes that of Bogomolny et 
al. [O]. These authors considered a single scatterer in a chaotic billiard and showed 
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Figure 2. Two-point correlation function R2{r) for different values of s = 2aA(e— £_). 
For comparison the GUE result is shown as well. 



that the fluctuations are chaotic. In the present contribution, we extend this study to 
arbitrarily many scatterers and also develop a completely different technique to derive 
the correlations. Generalizing the approach of Ref. ^H], the correlation functions are 
expressed in terms of a determinant. This determinant structure of the correlation 
functions is immediatly obvious in the supersymmetric formulation of the model due 
to the form of the Berezinian in eigenvalue angle coordinates. Moreover, an explicit 
Christoffel-Darboux formula is given for the kernel entering the determinant. 

By means of a saddlepoint appproximation, we rederive the density of states in the 
strong coupling limit and also find the real part of the one-point function. We show 
that the correlation functions in the bulk of the spectrum on the scale of the local mean 
level spacing are, for all couplings, of GUE type. 

Applying a WKB approximation to the kernel, the correlation functions are studied 
close to the band edges in the strong coupling limit, where the number of scatterers 
is large and the scattering potential is strong. The above mentioned saddlepoint 
approximation is not valid in this regime. Within the band the two-point correlation 
function shows a GUE behaviour, but approaching the band edges and proceeding 
towards the band tails the eigenvalues become more and more uncorrelated. This is 
exactly the fingerprint expected for a mobility edge and a localization-delocalization 
transition. We notice that a drastic depletion of the density of states this accompanies 
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this transition. Thus, the localization-delocalization transition found in the non-linear 
cr-model is of a different nature. In the latter, the average is over an ensemble 
of white-noise correlated impurities, while two averages are performed in the present 
model, one over the wavefunctions at the positions of the scatterers and another one over 
the billiard spectrum. The resulting models are therefore different. There is a kinetic 
term in the non-linear a-model and a diffusion constant in front of it. No analogy to this 
is present in the model discussed here, because the average over the billiard spectrum 
takes care of the kinetic term. 
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Appendix A. Properties of the Functions Generalizing the Hermite 
Polynomials 

It is useful to define the functions 

exp(-g 2 /2) 
<Pnl{£) = I Ar , u NL [e) 

^nl(s) = eXpi ~ £2/2) v NL (e) , (A.l) 



which reduce to the oscillator wave functions for L = or, equivalently, for a —>■ oo. 
We also introduce the operators 

A + = 4- - £ and A~ — + e , (A. 2) 

de de 

which act on the functions (|A.1|) according to 

A + (p NL (e) = - y/2(N+l)<p (N+1)L (e) 

A~ip NL (e) = + V2Nip {N -i )L (£) + —<pN(L+i)(e) 

A + i> NL {e) = - y/2(iV+ rj^ {N+1)L {e) - ^ N(L -i){e) 

A-^ NL {e) = + V2Ni/> {N - 1)L (e) . (A.3) 

These results extend the formulae for the oscillator wave functions by terms involving a 
change of the index L. We evaluate the action of the iterated operators A~A + and A + A~ 
using Eqs. (jA.3|) . properly combine terms and arrive at the second order differential 
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equations 

(J^- £2 + ( 2N + !)) Vnl{s) = ~yj2(N + l)<p (N+m+1) (e) 
( d 2 „ , „ A , , , L 



l J 2 e 2 + (2JV + 1)J = -^V2A^ M (L-i)(£) (A.4) 

These are no eigenvalue equations, because the functions on the left and the right hand 
sides have different indices. However, one can cast them into diffusion-type-of equations 
by introducing the fictitious time 

r = — In a such that a = exp(— r) . (A. 5) 

A straightforward calculation yields the equations 

" + {2N + 1} ) Vnl{£) = ~ 2 ^ nl{£) 
(J^2- £2 + ( 2N + !)) ^l(b) = +2^ NL (e) , (A.6) 
which involve the same indices on both sides. 
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